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It is shown that experimental values of the cross sections of inelastic low-energy neutron scattering 
on even-even nuclei together with the description of these cross sections in the framework of the 
coupled channel optical model may be considered as a reliable method for finding nuclei with a 
semimagic number (or numbers) of nucleons. Some examples of the application of this method are 
considered. 



I. INTRODUCTION 

During last two decades the problem of existence of 
so-called semimagic numbers of nucleons is widely dis- 
cussed (see e.g. [l|, where one can find a bibliography 
on this problem) . The existence of semimagic numbers is 
considered as a result of the appearance of a new gap in 
the single-nucleon level scheme by adding a pair of neu- 
trons (or protons) and its disappearance by adding an- 
other pair. Such an effect is believed to be caused by the 
interaction between valence protons and neutrons. For 
the first time the surmise about importance of the np- 
interaction was made by A.De Shalit and M.Goldhaber 
in the shell model in connection with the appear- 
ance of nuclear deformation. P.Federman and S.Pittel Q 
showed that the interaction between neutrons and pro- 
tons with strongly overlapped orbits may lead to the ap- 
pearance of nuclear deformation. At the same time they 
have shown that this interaction may essentially change 
the single-nucleon level scheme. This interaction becomes 
important when interacting neutron and proton have big 
radial or big and close orbital quantum numbers. There- 
fore semimagic quantum numbers are less steady than 
"classical" ones. 

As to methods of finding semimagic quantum numbers 
of nucleons it seems that shell model calculations taking 
into account the np-interaction of valence nucleons can 
not be regarded as reliable because the form and intensity 
of this interaction seem to be rather uncertain. There- 
fore semimagic numbers can be found by indirect ways, 
namely using a comparative analysis of some properties 
of neighboring nuclei, such as biding energies of nucle- 
ons, spectroscopic data (e.g. ^-factor values, change of 
the sign for the coefficient of E2-M1 mixture in electro- 
magnetic transitions 2 J — > 2+, 3^ — > 2^, 3^ — > 2^ of 
even-even nuclei etc). 

In the next section we show that the analysis of in- 
elastic scattering cross sections of low-energy neutrons 
on even-even nuclei gives a reliable method for finding 
semimagic numbers of nucleons. 

Section 3 presents some examples of the concrete ap- 
plication of this method. And section 4 contains some 
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conclusive remarks. 

II. NEUTRON INELASTIC SCATTERING 

Low-energy neutron data for even-even nuclei with 
A > 56 were successfully described [1, [f| by the au- 
thors of this paper in terms of the coupled channel op- 
tical model (CCOM). These data, taken for neutron en- 
ergies E n < 3MeV, included total and elastic scatter- 
ing cross sections, cross sections of inelastic scattering 
corresponding to the excitation of 2f level, angular dis- 
tributions for elastic and inelastic scattering. We used 
two-phonon (five-channel) approximation of CCOM for 
spherical nuclei, and three-channel rotational approxi- 
mation for nonspherical ones. Details of CCOM are de- 
scribed by E.S.Konobeevsky, I.V.Surkova et al [||. 

In our calculation we used nonspherical optical poten- 
tial with the real part of Woods-Saxon's form. The real 
part included the spin-orbital term and the symmetry 
potential. Radial dependance of the absorptive part was 
taken as the derivative of the real part of the form-factor 
since namely the nuclear surface is responsible for ab- 
sorption of low-energy neutrons. Geometrical parame- 
ters of the potential were the same for real and imagi- 
nary parts: the potential radius was fixed as R = roA 1 / 3 
with ro = 1.22 fm, the nuclear diffuseness parameter was 
initially chosen to be equal to 0.65 fm, but for some nu- 
clei (see below) it was somewhat changed. The spin-orbit 
interaction parameter V so was equal to 8 MeV. The real 
part of our potential included the isotopical term propor- 
tional to (N — Z)/A, so that the depth of the real part 
of the potential had the following form: 

N - Z 

V = V -V 1 —j-, (1) 

where V x = 22 MeV. The depths V (of real part) and W 
(of imaginary part) were two free parameters to fit. The 
values of the quadrupole deformation parameter /3 2 were 
taken from compillation [7] . 

The optimal description of all the totality of low-energy 
neutron data for even-even nuclei under consideration 
was achieved using practically the same values of model 
parameters Vq and W, namely, Vq = 52.5 ± 1.5 MeV and 
W = 2.5 ± 0.5 MeV. As to diffuseness parameter a, in or- 
der to get a good description of all the nuclei under con- 
sideration (including nonspherical ones) with the same 
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FIG. 1: Neutron inelastic cross sections for E n — 300 keV over 
threshold vs. the JVpjV n -product 



values of Vq and W, we had to change the diffuseness 
parameter for nonspherical nuclei from the initial value 
of 0.65 fm up to 0.70 — 0.75 fm. For the same reason we 
slightly diminished the value of a for some spherical nu- 
clei which were found to be magic or semimagic, taking 
a =0.55 - 0.60 fm for magic and semimagic nuclei and 
a ~ 0.50 fm for double-magic ones. 

These alterations in numerical values of a seem to be 
consistent (at least, qualitatively) with the experimen- 
tal data and theoretical view on the thickness of nuclear 
surface lay; on the other hand, they are important for 
finding semimagic numbers of nucleons (see below). 

Note, that using CCOM with theseparameter values 
we also obtained a good description 0, Q of the exper- 
imental data on s-, p- and d-neutron strength functions 
and potential scattering lengths for nuclei under consid- 
eration. 

The analysis of inelastic neutron cross sections led us 
to so-called iV p iV n -systematics 0,11]: we showed that the 
inelastic cross section with the 2^ level excitation, taken 
at the energy equal to 300 keV over the threshold and 
averaged over the energy range of 100 keV, seems to be 
a smooth function of N p N n — the product of valence 
proton and neutron numbers (or their holes) . The curve, 
presenting this function <Ti nc i(N p N n ) and shown with ex- 
perimental points in those publications, was obtained us- 
ing the least square method. The corrected and more 
precise version of this plot is given by fig.l of this paper. 
Note, that the curve presenting a function &i nc i(N p N n ) 
may also be well described as 

^incl 

= A + B^N P ~N~ n (2) 

with A =0.5 b and B =0.14b. Emphasize, that this ap- 
proximation is not good for big (> 120) values of N p N n . 

Since the iV p ^-product depends on the beginning 
and the end of the upper shells of the nucleus, the 
jVpiV n -systeniatics combined with CCOM-calculations 



may serve as a method of finding semimagic numbers of 
nucleons. The main scheme of this method is as follows. 

If the cross section of neutron inelastic scattering for an 
even-even nucleus is essentially less than for neighboring 
even-even nuclei (usually it is near by 0.5 b) and there- 
fore the corresponding experimental point on the plot of 
fig.l deviates from the curve <Ji nc i(N p N n ), we change the 
beginning (or the end) of valence shell in order "to re- 
store the agreement" . As a result of this operation, we 
obtain a new value of the valence nucleon product and 
thereby — a new number of nucleons corresponding to 
the beginning (or the end) of valence shell, i.e. a new 
magic number. This result may be considered as prelim- 
inary one; to confirm it we can see what is the result of 
the CCOM-description of this nucleus, namely, what is 
the value of diffuseness parameter: if the value of the dif- 
fuseness parameter is essentially less than 0.65 fm (0.50 
- 0.60 fm) it means that the preliminary result is con- 
firmed and a new semimagic number is found, if not — 
it is necessary to use another method to confirm or to 
reject the preliminary result. Such an important role of 
the diffuseness parameter in our method is based on the 
high sensitivity of the CCOM-calculation results to the 
diffusion parameter value. 

Thus, our method is expected to be efficient if neutron 
inelastic cross sections for the considered nucleus and for 
its "neighbors" are known with an appropriate accuracy. 
As we will show, in absence of such data it is also possible 
to obtain some results, but, may be, not so certain as in 
the main scheme of the method. 

Next section is devoted to practical applications of the 
method proposed here. 



III. SEMIMAGIC NUMBERS OF NUCLEONS 

Here, using the method described in the previous sec- 
tion, we prove the existence of neutron semimagic num- 
bers N s =38, 56, 64, and proton semimagic numbers Z s = 
40, 58, 64. We also consider possibilities of existence of 
some other semimagic numbers. 



A. N s =38 and Z s =40 

The case of N s = 38 seems to be the best illustration of 
our method application 0, Q (see also [§]). The neutron 
inelastic cross sections for neighboring isotopes 70 Ge and 
72 Ge differ by factor 1.7. So, following our scheme we 
have to assume that the JV p iV n -values for these isotopes 
are equal and 40 respectively. But it means that in 
70 Ge the neutron state pi/ 2 is shifted up, so that a con- 
siderable energy gap appears between it and low-lying 
state /5/2! this gap disappears by adding a pair of neu- 
trons, i.e. in 72 Ge. In other words, 7V S =38 is a semimagic 
number. Calculation of diffuseness parameter a for 70 Ge 
confirms this conclusion (the value a = 0, 65 presented in 
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4] does not take into account " the correction" of Vq and 
W; taking into account this effect gives a = 0.58). 

As to Z = 40 it seems that in distinction from N = 38 
the proton state P1/2 is not shifted up, and the energy 
gap appears between this state and the state g 9 / 2 (which 
is higher than p\i%). By adding a new pair of proton this 
gap disappears as a result of the interaction between va- 
lence protons and neutrons in states g 9 / 2 and 57/2- Thus, 
for protons we obtain Z s =i0. The standard (for our 
method) analysis of the low-energy neutron interaction 
with isotopes of Sr, Zr and Mo confirms this conclusion. 

Note, that the conclusion about N s and Z s follows 
immediately from the investigation by P.Federman and 
S.Pittel [3}. Later on many authors presented arguments 
in favor of a possibility of existence of semimagic numbers 
TV = 38 and Z = 40 (see e.g. flQj). But in distinction 
from those arguments our method gives evidence for the 
existence of these semimagic numbers. 



B. N s =56 and Z s =58 

To some extent the situation with semimagic numbers 
56 and 58 is similar to the case of numbers 38 and 40. 
The standard analysis of the experimental data on the in- 
teraction of low-energy neutrons with even-even isotopes 
of Sr, Zr and Mo 0] leads to the conclusion that for the 
isotopes with 56 neutrons the energy level of the neutron 
state 2d 5 / 2 is shifted down and the level of lg 7 / 2 -state is 
shifted up, so that a gap appears between these states and 
disappears by adding two more neutrons. In other words, 
N = 56 is a semimagic number. One of the first conclu- 
sion about a possibility of the existence of N s = 56 was 
made by V.A.Morozov d], who showed that at N = 56 
the coefficient of mixture of E2 and Ml in electromag- 
netic transitions 2^ —5-2^, 3^ — > 2^, 3^ — > 2^ changes 
the sign. And since such a change normally is treated 
as the presence of a filled subshell, he assumed that the 
existence of a semimagic number TV = 56 is quite pos- 
sible. Our result (being more definite) confirmed this 
assumption. 

In this connection we should mention recent works by 
Moscow University group [ll| in which it was in partic- 
ularly shown that the energy gap between neutron sub- 
shells in 96 Zr may achieve ~ 3MeV. As a result the au- 
thors assumed a possibility of existence of the semimagic 
neutron number N s = 56. This calculation may be re- 
garded as an additional confirmation for our result. 

As to a semimagic number of protons Z s =58 or 56 some 
authors (see e.g. [1]) assumed that such a number may 
exist. This assumption was confirmed in the framework 
of our approach. Trying to describe neutron data for Nd 
we discovered that neutron inelastic cross sections for 
some isotopes do not enter the 7VpAf n -systematics. For 
instance, from the point of view of this systematics the 
JVpiV n -value for the isotope 146 Nd must be not equal to 
40 (if its proton shell would start at Z—51), but 8 or 4, 
and these values correspond to Z s =58 or 56 (see fig. 2). 
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FIG. 2: Description of neutron date for Nd requires the 
iVpAT n -product to be equal not to 40, but to 4 or 8 which 
correspond to existence of semimagic numbers Z=58 or 56 

Considering isotopes of Ce (Z—58) and Ba (Z=56) 
shows that inelastic cross sections at the energy 300 keV 
above the threshold are somewhat lower for isotopes of 
Ce than for isotopes of Ba. For instance, for 140 Ce this 
cross section is equal to 0.60b but for 136 Ba — 0.66b. 
Diffuseness parameter a calculated by use of CCOM for 
isotopes of Ce is equal to 0.57 - 0.60 fm, but for isotopes 
of Ba — somewhat bigger (0.60 - 0.65 fm). 

Thus, our results confirm the assumption about ex- 
istence of the semimagic number Z s =58(56). At the 
same time they show that Z s =58 is more preferable than 
Z s =56. 



C. N S =Z S =64 

In distinction from the semimagic numbers considered 
above, the assumption about existence of the semimagic 
number 64 leads to the conclusion that some nuclei 
(namely 114 Sn and 146 Gd) must be double magic. 

In 1953 J.O.Rasmussen with colaborators ,12J payed 
attention to anomalies in energies of a-decay of nuclei 
with Z=64 at N « 82 and thus assumed that Z =64 is a 
semimagic number. Later on it was assumed [HI, US] that 
the subshell Z=64 exists in the region of rearearth nu- 
clei and disappears in consequence of the pn-interaction 
between fen /91- and fe 9 / 2 -states. Many authors (see e.g. 
[J, H, [H, EH, [IB]), considering the question about Z s , 
brought different arguments in favor of existence of such 
a semimagical number. In particular, the description of 
the neutron inelastic scattering for the isotopes of Sm in 
the framework of CCOM is possible only under assump- 
tion of the disappearance of the shell (1*77/2 and 2g? 5 / 2 ) 
gap at Z = 64 and AT > 80 @|. At the same time the 
value of diffuseness parameter a at the optimal descrip- 
tion of the neutron cross sections for 146 Gd is equal to 
0.50 fm, while for neighboring nuclei a=0.55 - 0.60 fm. 
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Such values of a give an evidence that 146 Gd is a double 
magic nucleus, i.e. Z—6A is a semimagic number. 

The situation with A s =64 is similar to that with 
Z s =64. Many authors 0, H GJ, [H| concluded that the 
existence of this semimagic number is quite possible. The 
analysis of neutron inelastic scattering on the isotopes of 
Cd in the framework of the iVpiVn-systematics confirmed 
that conclusion (just like such an analysis for Sm con- 
firmed the existence of Z S =6A). The value of the diffu- 
sion parameter for 114 Sn (a=0.50fm) gives an argument 
in favor of the conclusion that 114 Sn is a double magic 
nucleus (just like 146 Gd). 

Thus, we may conclude that the existence of semimagic 
numbers Z s = N s = 64 is firmly proved. 



D. Other possible semimagic numbers 

Besides the semimagic numbers of nucleons consid- 
ered above there are some "candidates" for joining the 
" semimagic community" . Unfortunately, the method de- 
scribed here cannot give so far any certain answer about 
their existence because in most cases the accuracy of nec- 
essary experimental data is insufficient. Nevertheless, the 
use of it can lead to some conclusions (at least, on the 
level of assumption). 

As the first example of such a " candidate" we consider 
Z s =76 [17]. Values of the neutron inelastic cross section 
for isotopes of Os (Z = 76) are essentially bigger than 
0,5 b and the diffuseness parameter values (corresponding 
to the best description of low-energy neutron data) are 
equal to 0.65. All that testifies against the existence of 
this semimagic number. 

It seems to be interesting to investigate the region of 
heavy [Z > 92) nuclei in connection with existence of 
semimagic nucleon numbers. But since for such nuclei 
the inelastic cross sections of low-energy neutron scat- 
tering are not known, it is impossible to use the main 
scheme of our method. However, one can try to draw 
some plausible assumptions (or at least some hints) from 
available data on inelastic neutron scattering on 232 Th. 
As a matter of fact there exist two sets of the experimen- 
tal data on neutron inelastic scattering for this nucleus. 
These cross section values corresponding to 2+ level exci- 
tation for the neutron energy 300 keV over the threshold 
(after averaging over the energy ran ge o f 100 keV) are 
equal to 1.13±0.21b and 1.22±0.22b Both of these 
values (corresponding to N p N n = 128) are essentially less 
than the value required by the A n A„-systematics. If we 
(ignoring the experimental errors) assume that the true 
value of the cross section considered here is situated be- 
tween 1.13 and 1.22 b, we obtain that the A p A n -product 
must be equal to 48 (instead of 128). Such a change may 
be correct only if some new semimagic numbers exist. 
We will try to find them. 

Generally speaking, the value of the A p A n -product 
equal to 48 may leads to 6 possible pairs of valence nu- 
cleon numbers N p and N n . But it is necessary that be- 
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The breach of one of these requirements would lead to 
drastic changes in the nucleon level schemes, in particu- 
lar, to the removal of classical magic numbers. (In princi- 
ple, such a reconstruction of the nucleon scheme is possi- 
ble. But there is no reason for it so far.). Therefore only 
three pairs (A p , A n )k (where k = 1, 2, 3), namely, (8,6), 
(6,8), (4,12) have to be taken into consideration. For 
each pair it is easy to get the nucleon numbers Z^ and 
Ak from which one has to count the values of relevant 
valence nucleon numbers: 



(4) 



Eqs.(|l]) contain 12 numbers Zy and Ak. Using these 
numbers one can form 12 pairs (Zk, Ak), and if both 
numbers of a pair are semimagic ones a necessary value 
of the ApA n -product (in our case it is 48) will be guaran- 
teed. It remains to prove that at least there is one such 
a pair among the pairs considered here. In this connec- 
tion note that one of the i?i-values is equal to 82 which 
is a classical magic number, so that in this case it is only 
necessary to prove that one of 136 or 148 is a neutron 
semimagic number. 

Unfortunately, it seems impossible to find a 
"semimagic pair" among those 12 ones for lack of nec- 
essary experimental data. Nevertheless, we shall try to 
draw some conclusions from available data (at least on 
the level of assumptions). Let us consider 12 nuclei con- 
sist of Z\ protons and Ak neutrons (Z^ and Ak belong 
to the same pair). 10 of these nuclei are not interest- 
ing for us because some of them are exotic ones (like 
228 Pb) and for others nothing is known about their in- 
teraction with low-energy neutrons. Remaining two nu- 
clei, namely gg 6 Cfi48 and g^Cmiso could be interesting 
from the point of view of our method. But we do not 
know anything about interaction of the isotopes of Cf 
with low-energy neutrons. So, only 246 Cm remains to 
be considered. For this and some other isotopes of Cm 
the neutron strength functions are known experimentally. 
The strength function So of s-neutrons for 246 Cm is equal 
to 0.45±0.15 while for neighboring isotopes 244 Cm and 
248 Cm S o =1.00±0.20 and 1.10±0.12 [l|. Such a big dif- 
ference between the So-values for Cm and its neighbor- 
ing even-even isotopes may be considered as an argument 
in favor of the assumption that A=150 is a semimagic 
number. The difference between the So-values for 246 Cm 
and 244 Pu is of the same order of magnitude, but we 
cannot say anything about the other " neighbor on the Z- 



line" 



, i.e. 



248 



Cf because its interaction with low-energy 
neutrons is not known. So, we can only say that such 
a situation may be considered not as an argument, but 
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only as a "hint in favor" of the assumption that Z=96 is 
also a semimagic number. 

Generally speaking, almost all the So-values in the 
(Pu-Cm)-region are equal to «1 The only exception 
is 246 Cm (So — 0.45 ± 0.15). Such a situation provokes 
us to assume that 246 Cm is a double magic (or double 
semimagic) nucleus. In other words, both Z = 96 and 
N = 150 are semimagic numbers. Emphasize, that this is 
only a statement about the possibility of existence of two 
new semimagic numbers and needs more rigorous proof. 
Unfortunately, the diffusion parameter value taken from 
our calculation of Sq cannot help because such calcula- 
tions are less sensitive to the diffusion parameter value 
than calculations of cross sections. 



IV. CONCLUSION 

The efficiency of the method for finding new semimagic 
numbers of nucleon, based on the analysis of low-energy 
neutron data, was demonstrated here by different ex- 
amples of its application. Using this method we ob- 
tained some already known results about existence of 
some semimagic numbers, confirmed assumptions about 
other semimags, for the first time proved the existence of 
the semimagic number Z s =58. Using scarce experimen- 
tal data we first pointed on the possibility of existence 
of semimagic Z — 96 and N = 150. Thus, this method 



seems to be rather efficient. 

It is necessary to emphasize that the most impor- 
tant part of the method described here is the N p N n - 
systematics, i.e., smooth dependence of neutron inelastic 
cross section on the iV p iV n -product. This dependence it- 
self is out of any doubt: there is no one serious contradic- 
tion with it. Therefore this method becomes less efficient 
in cases when one cannot address the -/V p 7V n -systematics, 
i.e. when neutron inelastic cross sections are not known 
or their accuracy is poor. 

Also we mention another essential element of our 
method having auxiliary character, namely the diffusion 
parameter value obtained from the requirement of opti- 
mal description of neutron data in framework of the cou- 
pled channel optical model. But this auxiliary element is 
efficient only if corresponding calculations are sensitive 
to this value. For instance, it works well for the cross 
section calculations, but does not work for calculations 
of strength functions (as it was mentioned above). Never- 
theless, even in such cases our method permits to obtain 
some results. 

Thus, the method described here may be consid- 
ered as one of the most efficient methods for find- 
ing new semimagic (or magic) number of nucleons. 
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